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\^ • Abstract 

o\ 

Topological defects must respect causality, a statement leading to restric- 

"-^ \ tive constraints on the power spectrum of the total cosmological perturba- 

["t I I tions they induce. Causality constraints have for long been known to require 

the presence of an under-density in the surrounding matter compensating 

the defect network on large scales. This so-called compensation can never 

be neglected and significantly complicates calculations in defect scenarios, eg. 

computing cosmic microwave background fluctuations. A quick and dirty way 

to implement the compensation are the so-called compensation fudge factors. 

CN| ' Here we derive the complete photon-baryon-CDM backreaction effects in de- 

^ ■ feet scenarios. The fudge factor comes out as an algebraic identity and so 

Q\ . we drop the negative qualifier "fudge" . The compensation scale is computed 

^ \ and physically interpreted. Secondary backreaction effects exist, and neglect- 

Qh| ing them constitutes the well-defined approximation scheme within which one 

O ' should consider compensation factor calculations. We quantitatively assess 

-*— > ' the accuracy of this approximation, and conclude that the considerable pains 

d • associated with improving on it are often a waste of effort. 



o 
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I. INTRODUCTION 

When a defect network is formed, causality and energy conservation demand that there 
must be a compensating under-density in the background matter and radiation. This com- 
pensation is exactly anti-correlated with the defects and is of a comparable intensity. It acts 
as a source for the gravitational potential which in turn drives radiation perturbations. For 
this reason the compensation cannot be ignored in CMB calculations. In previous work [jl],|^ 
the compensation has been included by use of compensation fudge factors. These ensure 
that the overall perturbations have a large-scale behaviour consistent with the causality 
constraints 0,^. However their exact form was never justified by an analytical identity, and 
hence the words "fudge factor" qualifying them. 

In this paper we show that by choosing a suitable gauge it is possible to derive analytically 
an expression for the compensation which contains a term in the form of the previously 
discussed fudge factor. Along with the compensation factor we find that there are also 
terms representing the backreaction from baryons, CDM and radiation which cannot be 
predicted a priori. We find however that these terms are sub-dominant compared to the 
defects for any reasonable values of the Hubble constant and the baryon fraction today. 
The physical implication is that it is the defect rather than any other component which 
dominates the spectrum of radiation perturbations. Hence by choosing a suitable defect it 
is possible to have a range of spectra which do not necessarily have the characteristic out of 
phase spectra associated with defect calculations, as shown in [Q. 

This paper is organised as follows. In Section H review the equations for a system of 
photons, baryons, CDM, and defects, in the tight-coupled approximation. We present a trick 
for easily considering a defect component within a multi-fluid formalism, such as the one 
presented in Kodama and Sasaki ||^. We show how the various gauge-invariant formulations 
correspond to nothing but different choices of full gauge fixing. We look into all popular 
choices of gauge. We argue in favour of the flat slicing gauge for the discussion of feedback 
mechanisms, the compensation, or the causality constraint. Then in Section |I| we present 
an algebraic manipulation which allows splitting compensation into two terms, one of which 
is purely a compensation fudge factor. We evaluate then, in Section |IV|, the quantitative 
impact of dropping the terms other than the compensation factor term. We find that one 
would need rather extremely values of the Hubble constant and baryon fraction for these 
extra terms to have much qualitative effect. In the concluding section we digress on the 
metaphysical implications of this results, and the practical application which we intend to 
give it. 

II. THE BASIC EQUATIONS 

We consider the epoch before recombination when it is a good approximation to treat all 
the cosmic components as a fluid. This is the so-called tight-coupled approximation, which 
we shall use to develop all our arguments. We will consider a scenario where the Universe is 
made up of radiation (corresponding variables labelled by 7), baryons (6), cold dark matter 
(c), and a defect component (s) where the baryons and radiation are tightly coupled (see 
eg. M for a quantitative definition). 



Whenever we intend to show the generahty of our results we shall also consider a set 
of extra generic components, denoted by a, which may be neutrinos or whatever personal 
taste requires. In order to extract intuition from our arguments we shall sometimes assume 
that only defects and radiation are relevant, a statement valid deep in the radiation dom- 
inated epoch. The actual calculations will however always be valid considering the other 
components. This is necessary for generality, since matter-radiation equality may (and in 
fact usually does) happen before recombination. 

We use the gauge-invariant formalism in all the guises discussed in Kodama and Sasaki 
1^ (KS from now on). These different formulations, one should point out, correspond to 
different choices of gauge. A set of a priori gauge dependent variables defined in a fully fixed 
gauge is of course gauge invariant. The different possible gauge-invariant contrast variables 
correspond to nothing but the density contrast as measured in different fully fixed gauges. 
There is sometimes an inappropriate feeling that "gauge-invariant" and "gauge-dependent" 
methods are two separate tool boxes. They are in fact one and the same thing. The only 
exception to this statement is the synchronous gauge 0. This choice of gauge leaves a 
residual gauge freedom and therefore variables defined in synchronous gauge can never be 
related algebraically to gauge invariant variables (see for a good dictionary). 

In this paper we adopt a multi-lingual attitude towards cosmic perturbation formalisms. 
We shall use a variety of density contrast variables A^ for a generic component a. These will 
be indexed in the following way. No extra index after the component index a denotes density 
contrast in the a-component rest frame (not in the total rest frame which we found a mess 
in the presence of defects). An index s denotes the density contrast in the Newtonian slicing 
(where the perturbed cosmological flow appears to have no shear). An index g denotes the 
density contrast in the fiat slicing (where the equal time slices have no scalar curvature). 

We consider all these different gauges in order to connect with previous work. However 
the main remark in this and the next section is that by choosing the fiat slicing contrast 
variables two desirable (at least for defect practitioners) features may be achieved. Firstly, 
one gets rid of potential time derivatives in the equations for the radiation (and also for 
baryons and CDM, but this can also be achieved in other gauges). This is a major tech- 
nical improvement on the formalism. It allows performing all calculations invoking only 
defect stress energy components and not their time derivatives. Structure functions for time 
derivatives are notoriously noisier in defect simulations. 

Secondly, as shown in the next Section, in the fiat slicing gauge the radiation backre- 
action naturally separates into two terms. The first is required by the Traschen integral 
constraints [Q and cannot be set to zero in any approximation scheme, otherwise causal- 
ity is grossly violated. We will however find an exact expression for this term made up of 
a factor (independent of the perturbation variables) times a set of defect variables. This 
factor happens to have the same form as the "compensation fudge factor". Furthermore 
the "compensation scale", left undetermined in compensation factor calculations, can be 
computed. Hence fudging is an exact equality made obvious in the fiat-slice gauge. The 
second backreaction term is truly unpredictable, but we will be able to show that it is not 
required by causality, and it is qualitatively sub-dominant. 



A. A defect component in a multi-fluid formalism 

We shall add to the multi-fluid formulation of Kodama and Sasaki a defect component. 
This can be best implemented by noticing that defects have no background stress-energy. 
We may then regard defects as a fluid for which the background energy and pressure are zero, 
the perturbation variables are infinite, but the product of the background and perturbation 
variables is finite and equals the defect stress-energy. 

More mathematically let the defects stress-energy tensor B^j, be a pure scalar so that it 
may be written as 



eoo = p' 




Go^ = hv' 




Qij = p'Sij + {kikj - 


- ls^,k')U■ 



(1) 

Then let us consider a component a = d with pd = Pd = 0, infinite perturbation variables 
(eg. 6d = C)o), but finite products of the two. From the way perturbation variables are 
defined in KS from the stress-energy tensor we can then write 

a'^pdSd = p' 
a^{Pd + Pd)vd = kv^ 

a^Pdifd = ^'n^ 

aVdHi = p' (2) 

Because the background stress-energy of defects is zero, defect variables are gauge-invariant 
by themselves. However care must be taken when identifying KS defect variables with defect 
variables. For instance, it may happen that a gauge invariant density contrast variable is 
given by a combination of defect variables. Using Eqns. ^ we can find the identifications: 

aVdArf = p' + 3hv' (3) 

a^pd^sd = a^Pd^gd = p' (4) 

For all other variables there is no ambiguity, as the extra terms required to turn gauge- 
dependent variables into gauge independent ones simply vanish. For instance 

a^iPd+PdWd = o?{pd + Pd){vd - Hr/k) = kv'' (5) 

a2p,n^=p^ (6) 

a^p.nj = k^W (7) 

The conservation equations for the defect component may be written as; 

p' + /i(3/ + p') + k^v' = (8) 

v' + 2hv' -p' + h^W = . (9) 

o 

which can be derived from the conservation equations in KS with the identifications made 
above. The gauge-invariant potentials $ and \1/ can be obtained from the Einstein's equations 



in KS. These are now sourced by a total density contrast and anisotropic stress which 
contains defects. We choose to separate the defects from all other components. Hence in all 
formulae in KS containing totals the following replacements should be introduced 

a^Ar -^ a^P^T + p' + ^hv' = a^Y. P"^" + P' + '^hv' (10) 

a^{p + p)Vt ^a^{p + p)Vt + kv' = c? ^(p„ + p^)V^ + kv' (11) 

a^n? ^ a^n? + k^W = a^Yl P^^l + '^'n'^ (12) 

Bearing this in mind, the Einstein equations in the presence of defects may now be read off 
from KS as 

A;2$ = 47r(aVAT + p' + 3hv') (13) 

^ + ^ = -8n(a'^ + U^) (14) 

In the scenario we are considering the total density contrast, putting defects aside, is given 
by 

pAt = PbAb + p^A^ + pcAc (15) 

The fluids viscosity 11^ is entirely due to the photons brightness quadrupole |Q]: 

n = y 92 (16) 

and can be set to zero in the tight-coupling limit. 



B. The Newtonian slicing equations 

The Newtonian slicing equations for the radiation are what leads to the Hu and Sug- 
yama (HS) formalism pj. During tight-coupling the photon system is described in HS by 
the monopole and dipole components of the brightness function, Bq and 0i. In the fluid 
description this corresponds to the Newtonian slicing variables 

eo = A,^/4 (17) 

Oi = V^ (18) 

It can be checked that, with this identification, the conservation equations in KS for radiation 
become the HS equations: 

where R = |— is the scale factor normalised to 3/4 at photon-baryon equality. It is for 
these equations that HS propose a WKB solution, which was used in the study of defect 



Doppler peaks in [0] . If one uses this gauge for the radiation one must however change to 
the comoving gauge before computing the potentials \l/ and $. This can be easily done by 
means of 

A^ = 4(^00 + /i^) (20) 



C. The comoving gauge 

In Hu and White (HW) |^ the issue of backreaction is addressed in the comoving gauge 
(the word gauge being replaced by "representation"). A temperature variable is defined such 
that 

r = A^/4 = eo + h^ (21) 

and an horrible set of equations is derived for them. The comoving gauge has the advantage 
that it is the natural gauge for representing the baryons, since in tight coupling baryons and 
photons share the same rest frame. As shown in HS and HW, the baryons' density contrast 
Ah and velocity VJ, satisfy the conditions 

A, = ^A, (22) 

Vt = V^ (23) 

which can be rewritten by defining the entropy as s = A;, — (3/4) A^, and rewriting the first 
equation as s = 0. 

The comoving gauge is also the gauge where the Traschen integral constraints are written 
1^,^. In an expanding Universe a set of energy conservation laws apply to perturbation 
variables [Q. When applied to a Universe which is initially unperturbed, and then causally 
made inhomogeneous, these constraints translate into a stringent requirement on the large 
scale power spectrum of these perturbations |^. This requirement is roughly that the power 
spectrum of the total energy perturbation in the comoving gauge goes like k'^ for small k. 
In the presence of defects the energy density subject to this law is: 

U = a'^pAr + p' + 3hv' (24) 

which is also the source of the gauge-invariant potential $. Hence one can rephrase the 
causal constraint as the requirement that $ goes to a constant at low k (white-noise). 

The lack of superhorizon correlations in the defect network requires the power spectrum 
P{p^) to have a white noise low k tail. Energy conservation Eqn. (^ requires that f^ 
also have a white noise low k tail. Hence the causal constraint entails the need for the 
compensation: a low k white-noise tail in the power spectrum of non-defect matter. The 
compensation must be exactly anticorrelated with the defects' tail. The quantity to be 
cancelled is p* -|- Shv"^, and not just p'^. The density forced to have a k'^ power spectrum is 
U and not just a combination of At and p'*. 



D. The flat-slicing gauge 



We can also define a temperature perturbation in the fiat-slicing gauge |TB]: 

Ao = Oo + $ = A,^/4 (25) 

In this gauge the photon equations are: 

As announced before one gets rid of the potential time derivatives in this gauge. We will 
also find it useful to represent CDM in this gauge, so that CDM equations are: 

V^ = -hVc + fc^ (27) 

There is a good mathematical reason why this gauge may be better for discussing causality 
and compensation issues. This is the gauge where the fluid equations of motion more 
resemble Minkowski space-time equations of motion. Hence energy variables in this gauge are 
akin to the pseudo-energy usually defined in the synchronous gauge, and used to introduce 
the compensation |[T^ , p!5[| . 



III. EXACT COMPENSATION FACTORS 

We start with an algebraic remark. The source for the Einsteins equations are energy 
density variables in the comoving gauge. In the scenario we are considering 

A;2$ = 4:7i{a^p{nbAb + n^A^ + ficAc + fi^A^) + p" + 3hv'} (28) 

$ + ^ = -Svrff (29) 

where a represents any other component we may have forgotten, with a sum over a implied, 
if need be. Now let us express baryons in terms of photons by means of Eq. ( p2D written as 



3 
Ah = -A^ + s (30) 



and write all other sources in the fiat-slice gauge: 



A, = 4(Ao-$ + /i^) (31) 

A, = A,, + 3(^/1^-$) (32) 

A, = A,, + 3(1 + Wo.) (h^ - $) (33) 



With these rearrangements, the first Einstein equation becomes 

a(Ac5 - 3$ + 3^) + n^l^Aag + 3(1 + «;„)(-$ + ^))) + p^ + 3hv') (34) 

The source term can now be spht into 3 components: 

k^<l> = S + Si + S2 (35) 

where 

S = 47T{a^pbS + p' + 3hv') (36) 

S^ = -47raV(^7(l + R) + 3fic + 3(1 + Wa)^^)* (37) 

^2 = 47raV(4fi^(l + R){Ao + hOi/k) + Vl^iAg^ + 3hVjk) + ^^(Ag^ + ShV^/k) (38) 

S is made up of sources which drive the radiation-baryon-CDM system but which are external 
to them. We call it the external source. This may be a topological defect. An entropy 
perturbation may also be regarded as external since it evolves independently of all other 
perturbations (according to s = 0). There is a fundamental difference between defects and 
entropy perturbations. Entropy perturbations satisfy s = 0. Defect sources, on the contrary, 
satisfy Eqns. (|) and (|^). 

The photon-baryon-CDM system is also driven by a backreacting term, here split as 
Si + S2- This reflects the fact that baryons, photons and CDM are driven by a potential 
which they are a source of. There is therefore a (linear) feedback effect which jeopardises 
for instance the use of the WKB solution in HS for defects. One could hope that defects are 
the main driving force, and try to neglect backreaction. However this back-reacting term 
incorporates the compensation. Setting Si + S2 to zero is therefore an approximation which 
can never make sense, as it would imply a gross violation of the causality constraint. The 
potential $ power spectrum would diverge like l//c^ at small k rather than go to white noise. 

However, in the fiat-slice gauge we have an algebraic bootstrap which one may hope 
already reflects most of the physical feedback mechanism. This bootstrap is created by 
the term Si. Let us first set 6*2 = 0. Then all the backreaction is predictable and fully 
determined by the defect sources. 5*1 may be passed to the left hand side of the Einstein 
equation (^) and be incorporated in an equation where the external sources are simply 
multiplied by a factor independent of photon, baryon or CDM variables. More precisely 



A;2$ = 47r7c(aVbS + p' + 3hv') (39) 



where 



xl = lihr])\4n^{l + R) + 3fie + 3(1 + W7„)fi,) (41) 

The backreaction encoded in Si is not an independent feedback mechanism operating in 
the fluid and imprinting a fixed signature in the photons' power spectrum for any defect 
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theory. From this term we can never expect to derive an out-of-phase signature as the one 
attributed to defects in |^. The backreaction contained in this term is fully driven by the 
defects alone, and can be made to behave in whatever way we want by properly designing 
the defect. It is not surprising that by considering only this backreaction effect we can place 
the primary Doppler peak anywhere, including the adiabatic and out-of-phase positions . 
On the other hand by considering this term one is already taking into account the causality 
constraint. The potential $ according to the new equation ( |5^ already goes to white noise 
at small k. 

If 5*2 7^ then one must add to equation (|5D|) an extra source term, so that 



e^ = 7c(5 + ^2) (42) 

This is an exact expression. The compensation factor has appeared as a result of algebra 
and the compensation scale is a well defined quantity dependent only on the unperturbed 
cosmological expansion dynamics. The compensation factor approximation is now the claim 
that we can set to zero 82- This is a well defined statement which we should be able to assess 
quantitatively. The term in S2 is the truly unpredictable backreaction. The 5*1 + 5*2 split 
has allowed us to separate what is truly a problem and what is not. By doing so we have 
implemented an approximation scheme (52 ~ 0) where we may avoid the feedback problem 
without immediately doing something stupid, like violating the causality constraint. 

Physically what this algebraic manipulation amounts to is the realization, made obvious 
in the flat-slice gauge, that the compensation is made up of 3 terms. One is the energy 
perturbation as it appears in a gauge where equal time surfaces appear to have no curvature. 
The other two are a potential perturbation describing the curvature of these slices, and a 
velocity term. The potential term is caused by the defects as well, and by dropping all other 
contributions we obtain a non pathological approximation scheme where the compensation 
is fully gravitational and perfectly correlated to the defect network. 

The compensation scale Xc varies from a/6 in the radiation epoch to vTS in the matter 
epoch. The compensation scale depends purely on the expansion kinematics and is affected 
by the matter radiation transition. It can be written as 

Xl = 3r]\h'^ -h) = Una^p + p)rf (43) 

and ii h = a/rj then x^ = 3Q;(a + 1). 

IV. QUANTITATIVE ARGUMENT IN FAVOUR OF COMPENSATION 

FACTORS 

The question remains of how good an approximation setting 6*2 = is. We address this 



question quantitatively. For definiteness we use the source defined in |10 



^ _ 1 sinAA:?7 
P - ^1/2 ^^^ l^^J 

ff = (45) 

bearing in mind that this ansatz may preclude arbitrary shifts in peaks' positions. This 
property is far from general, as shown in IlllO]- This issue is beyond the scope of this 
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paper, but will be addressed in a future publication |]T^. We then solve equations ( p6D for 



radiation, (27) for CDM, (H) and (P) for the remaining defect variables, and (13) and (|T4| 
for the potentials (with ( [T3| ) rewritten as in (|^)). The baryons are solved implicitly by the 
tight-coupled conditions (|2^) . By including 5*2 as a source of the potential $ in (^) one 
is considering the full backreaction effects, due to baryons, radiation, and CDM. One thus 
obtains the exact solution to this problem. By dropping the term in Age and Vc in 5*2 one 
neglects the effects of CDM fluctuations on the CMB fluctuations. By setting 5*2 = one 
neglects the effects of backreaction altogether. In the last approximation the source is truly 
external, and is compensated purely by defect gravitational effects. This is the compensation 
factor approximation. 

We have solved this problem for various values of the Hubble constant and baryon content 
of the Universe. These are parameterised by h, so that the Hubble constant nowadays is 
Hq = 100/iKms~^Mpc~^, and ^h = Pw/Po, the baryons density fraction nowadays. We 
expressed our results in terms of the effective temperature Gq + \I/ = Aq + (^ — $) and 
the dipole 6i at last scattering t] = r]^. This is because these are the quantities which are 
then projected onto Ci's by free-streaming after the last scattering surface 0. We have 
solved the full problem and compared the full answer with the effect of dropping CDM, and 
dropping all CDM-photon-baryon backreaction. In all plots these 3 calculation schemes will 
be represented by dash, dots, and lines, respectively. 

The results for a source with A = 1 are plotted in Figure |^, and we now comment on 
them. We have chosen extreme values for Q^ and h in order to emphasise the point we wish 
to make. For popular values of Qb and h the compensation factor approximation works very 
well. Clearly one needs rather high values for both h and Qb for the compensation factor 
approximation to become gross. 

As h increases the time between equality and last scattering increases. As a result CDM 
fluctuations have time to start to grow while they still can interfere with the tightly coupled 
radiation. As a result the CDM contribution cannot really be neglected in scenarios with a 
large Hubble constant. 

As fib increases the so-called acoustic signature may be imparted on the peaks. This 
is an asymmetry in amplitude between odd and even peaks resulting from baryons shifting 
the zero level of the oscillations. After squaring, the peaks will appear alternately big and 
small. It is interesting to note two things. First the acoustic signature appears already in 
the compensation factor approximation, although less pronounced. Secondly, if one is to pay 
attention to detail, then CDM is as important as the baryons in imprinting the full acoustic 
signature. 

Another feature present in the spectra is the shoulder preceding the peaks, ubiquitous in 
defect scenario. This is not present in Aq and is due to the gravitational redshift term \1/ — $. 
At large scales Aq ~ and then goes negative. The potential term \l/ — $ is white noise 
and positive as /c — >• 0, then goes to zero. This induces a pre-peak which is not acoustic, but 
merely a gravitational redshift effect at last scattering. 

All in all backreaction seems to be negligible in the qualitative study of defects. His- 
torically it has been known that backreaction in defect theories can never be neglected, as 
a result of the causality constraints. However we have now shown that backreaction can 
be exactly split into two terms. One is purely gravitational, makes sure that the causal- 
ity constraints are satisfied, and can be predicted a priori from the defects by means of 
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a compensation factor. The other cannot be predicted a priori but it is not required by 
causahty. It merely reflects the baryons, photons, and CDM trying to make a nuisance of 
themselves, rebelling against the driving force of the defects, trying to imprint a signature 
they are allowed to imprint whenever no driving force is present. However for all reasonable 
values of Qb and h this secondary backreaction is quantitatively sub-dominant. When driven 
by defects the radiation/baryon/CDM feedback effect is normally weaker than the external 
force they are subject to. Therefore, it turns out that in defect scenarios, once compensation 
factors are taken into account, backreaction is precisely something which can be neglected, 
certainly in any qualitative discussion of defect perturbations. 

V. WASTE NOT - A DEFECT FACTORY 

In this paper we stressed the technical difficulties induced when considering compensation 
in defect calculations but have shown how they can be effectively bypassed by means of 
compensation factors. The compensation results from a feedback mechanism which operates 
in the photon, CDM, baryon system when an external driving force is applied to them. Some 
of this backreaction is unpredictable in the sense that by specifying the driving force one 
does not specify the backreaction before the system of differential equations describing the 
whole system is evolved. However the dominant compensation term is always predictable 
in the sense that all one needs to do is to multiply some appropriate combination of defect 
stress energy components by a compensation factor. 

The practical implication of this result is of course not to simplify the numerics of evolving 
the tight coupled equations; this problem is straightforward enough for simplification to be 
necessary. However there is another, more metaphysical, side to compensation considerations 
in the literature: this is the belief that backreaction effects force the Doppler peaks to be 
out of phase. This could only be a theorem if the compensation was the result of a self- 
regulatory feedback mechanism operating in the fiuid, regardless of the details of the force 
driving it. The compensation described by the compensation factor, which we showed to be 
dominant, is precisely the opposite of this. It represents a feedback effect which is directly 
connected to the driving force, and which has a scale which can be manipulated by tuning 
the scale of the driving force. It was in fact shown in [^ that in the compensation factor 
approximation, a causal source could be designed so as to shift the Doppler peaks about, 
placing them on the adiabatic position or anywhere to its right. 

The fact that the compensation factor approximation also seems to reproduce other 
effects, such as the acoustic signature, then provides us with the practical application of 
this work. If the compensation factor approximation works one can write down a solution 
expressing the monopole and dipole in the photons at last scattering as a function of a 
defect structure function. We can then invert this expression and instead write down a 
defect structure function which could give us a given monopole and dipole at last scattering. 
This trick can then be converted into a defect factory, allowing intelligent guesses of causal 
sources which exhibit effects under study. In particular one may use this defect factory to 
produce confusing defects: causal sources exhibiting allegedly inflationary signatures, such 
as the acoustic signature. Naturally once an educated guess has been made, one should then 
go and evolve the full system of tightly coupled equations, or even better solve a Boltzmann 
code for the proposed source. However there is no reason why one should go the complicated 
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FIG. 1. We plot k^'^{@Q + ^) (curves which go negative first) and /c^'^Bi with the fonowing 
conventions: the compensation factor approximation in lines, photon-baryon backreaction included 
in dots, all backreaction included in dash. From left to right Jlfe = 0.005, 0.05, 0.5. From bottom 
to top h = 0.1,0.5,0.8. 
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way when guessing sources. In a future publication we shall use this approach in the study 



of confusing defects Hi? 
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